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Abstract 

The celebrated Zeilberger algorithm which hnds holonomic recurrence equations for 
definite sums of hypergeometric terms F(n, k) is extended to certain nonhypergeometric 
terms. An expression F(n, k) is called hypergeometric term if both F(n + 1, k)/ F(n, k) 
and F(n, k + 1)/F(n, k) are rational functions. Typical examples are ratios of products of 
exponentials, factorials, T function terms, binomial coefficients, and Pochhammer symbols 
that are integer-linear with respect to n and k in their arguments. 

We consider the more general case of ratios of products of exponentials, factorials, 
r function terms, binomial coefficients, and Pochhammer symbols that are rational- linear 
with respect to n and k in their arguments, and present an extended version of Zcilbcrger's 
algorithm for this case, using an extended version of Gosper's algorithm for indefinite 
summation. 

In a similar way the Wilf-Zeilberger method of rational function certification of integer- 
linear hypergeometric identities is extended to rational-linear hypergeometric identities. 

The given algorithms on definite summation apply to many cases in the literature 
to which neither the Zeilberger approach nor the Wilf-Zeilberger method is applicable. 
Examples of this type are given by theorems of Watson and Whipple, and a large list of 
identities ( "Strange evaluations of hypergeometric series" ) that were studied by Gessel and 
Stanton. It turns out that with our extended algorithms practically all hypergeometric 
identities in the literature can be verified. 

Finally we show how the algorithms can be used to generate new identities. 

Reduce and Maple implementations of the given algorithms can be obtained from 
the author, many results of which are presented in the paper. 



1 Hypergeometric identities 



In this paper we deal with hypergeometric identities. As usual, the notation of the generalized 
hypergeometric function p F q defined by 
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is used, (a)k = T< Y(k)' > denoting the Pochhammer symbol or shifted factorial. The numbers a k 
are called the upper, and the lower parameters of p F q . 

The coefficients A^ of the generalized hypergeometric function have the rational term ratio 
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If an expression Ak has a rational term ratio Ak+i/Ak, we call Ak a hypergeometric term 
or closed form. Note that any hypergeometric term essentially has a representation as the 
ratio of shifted factorials (over (D), and its generating function is connected with a generalized 
hypergeometric series. 

The classical reference concerning generalized hypergeometric series is the book of Bailey 
(1935) containing a huge amount of relations between hypergeometric series some of which 
represent the value of certain hypergeometric functions at a special point (mostly x = 1 or 
x = — 1) by a single hypergeometric term. We will be concerned with this type of identities, 
and Table |] is a complete list of all such hypergeometric identities found in Bailey's book. 
Here n € IN is assumed to represent a positive integer so that the hypergeometric series 
with upper parameter — n are terminating. All other parameters involved represent arbitrary 
complex variables such that none of the lower parameters corresponds to a negative integer. 
With a method due to Wilf and Zeilberger, and with an algorithm of Zeilberger, many of 
these hypergeometric identities can be checked. It turns out, however, that for some of these 
identities both methods fail. We give extensions of both the Wilf-Zeilberger approach, and 
the (fast) Zeilberger algorithm with which all above identities can be handled as well as a 
large list of identities that were studied by Gessel and Stanton (1982). 
Our extensions therefore unify the verification of hypergeometric identities. 



2 The Gosper Algorithm 

In this section we recall the celebrated Gosper algorithm (Gosper, 1978), see also Graham, 
Knuth and Patashnik (1994). 

The Gosper algorithm deals with the question to find an antidifference Sk for given ak, i. e. a 
sequence Sk for which 

afc = s k - Sfc-i (2) 
in the particular case that Sk is a hypergeometric term, therefore 

— — is a rational function with respect to k , (3) 
i. e. Sk/sk-i £ Q(&). We call this indefinite summation. 

Note that if a hypergeometric term antidifference Sk exists, we call the input function 
Gosper-summable which then itself is a hypergeometric term since by (^) and (||) 
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Ofc-l Sfc_l — Sk-2 1 — Vk 

is rational, i. e. Uk,Vk £ Q[fc] are polynomials. 

Now, Gosper uses a representation lemma for rational functions to express akja-k-i in terms 
of polynomials. 

The idea behind this step comes from the following observation: If we calculate a& from 
s k = (2k)\/k\, e.g., we get 

ak = st _ ^ = (f _ M = mm - 1) - *, . <»z*>! = m - 3) . " 

i. e. a product of a polynomial pk = k(4k — 3) and a factorial term bk = jEf~^ f° r which 
bk/bk-i = qk/ r k is rational, and therefore and can be assumed to be polynomials. 
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Table 1: Bailey's hypergeometric database 
page Theorem Identity 



2-3 



Vandcrmondc 



Gaufi 
9 Saalschiitz 

9 Kummer 

11 GauB 

11 Bailey 

13 Dixon 



16 



Watson 



Whipple 
16 Whipple 
DougalPs 



26 



Theorem 



25/27 Dougall 



28 Whipple 

30 Bailey 

30 Bailey 

30 Bailey 

30 Bailey 



3-^2 

2-^1 
3F2 



3F2 



3-^2 



7F 6 



5F4 
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3-^2 



4^3 



4-F3 



a , b 

c 



1 = 
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a , 6 , — n 
c,l + a + fe — c — n 



T(c-a)T(c-b) 

(c - a) n (c - b) n 



a , 6 
l + a-b 



-1 = 



y (c) n (c-a-6) n 

(i + o)_ 6 r(i + o-6)r(i + o/2) 



(l + a/2)_ 6 r(i + o)r(l + o/2 - 6) 
r(l/2)r((a + 6+l)/2) 



a ,6 1 \ 

(a + 6 + l)/2 2 J ~ r((a + l)/2)r((6 + l)/2) 

A r( c /2)r(( c + i)/2) 



a , 1 — a 
c 



a , 6 , c 
1 + a — 6, 1 + a — c 



2/ r((a + c)/2)r((l-a + c)/2) 

(l + o)_ c (l + o/2-6). 
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(l + a/2)_ c (l+a-o)_ c 
r(l + a/2)r(l + a - b)r(l + a - c)r(l + a/2 - b - c) 
r(l + a)r(l + a/2 - 6)r(l + a/2 - c)T(l + a - b - c) 

r(i)r(i^)r(i+§+^)r( 1 - a - 2 b + 2c ) 



a , b , c 
(a + 6+ l)/2 ,2c 



a, 1 — a, c 
, 1 + 2c - e 



1 = 



; r(i^)r(i^)r(±^^)r(M±2c) 

7r2 1 - 2c r(e)r(l + 2c-e) 

r(^)r( a+1 + 2c - e )r(i^f^)r( 2 + 2c 2 ' a - e ) 



a , l + a/2 , 6 , c , d , 1 + 2a — b — c — d + ra , — n 
a/2 , 1+a— , 1+a — c,l + a — d , b+c+d—a—n , 1 + a + n 

(1 + a) n (1 + a - b - c) n (1 + a - b - d) n (1 + a - c - d) n 
(1 + a - 6)„(1 + a - c)„ (1 + a - d) n (1 + a - b - c - d) n 

a,l+a/2,c,d,e \ _ (l + a)_ e (l + a-c-d)_ e _ 

a/2, 1+a-c, 1+a—d, 1+a — e J (l+a-c)- e (l+a-d)- e 
r(l + a - c)r(l + a - d)r(l + a - e)r(l + a- c- d-e) 
T(l + a)r(l + a-d- e)r(l + a - c - e)r(l + a - c - d) 

(l+a)_ e r(l+a-d)r(l+a-e) 



a, l+a/2, d, e 
a/2, 1 + a — d, 1 + a — e 



a , 1 + a/2 , — n 
a/2 , w 
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a , 6 , — n 
l + a- 6,l + 26-n 



J \l+a-d)- e r(l + a)r(l + a-d-e) 
(w — a — 1 — n)(w — a) n _i 

(w)n 

(a-26)„(l + a/2-6)„(-6)„ 
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a , 1 + a/2 , b , — n 
a/2 ,1 + a- b , 1 + 26- n 



a, l + a/2, 6, — n 
/2, l + a-6,2 + 2o-ri 
3 



1 = 



(l + a-6)„(a/2-6)„(-26)„ 

\ = (a-2b) n (-b) n 
J (l + o-6)„(-26)„ 

(a-26-l)»(l/2 + a/2-6)„(-b-l)» 
(l + a-6)„(a/2-6-l/2)„(-26-l)„ 



Gosper shows then that such a representation with the property 



gcd (qk, ^k+j ) = 1 for all j G No (4) 



generally can be found and gives an algorithm to generate it. 
Therefore we have for ak the relation 



a-k Pk Qk 

(5) 



a-k-i Pk-i r k 

Pk corresponding to the polynomial and (<Zfc, r fc) to the factorial part of a^. 
Gosper finally defines the function by the equation 



Qk + l , /R \ 

Sk = Jk a k (6) 

Pk 



for which one sees immediately that 



n Pk S k Pk S k Pk 

Jl 



Sk 

Qk+l Ofc Qk+l Sk — Sk-i qk+i — - 



Sk-l 

is rational. Using @), Gosper proves the essential fact that fk is a polynomial. 
It follows from its defining equation that the polynomial fk satisfies 

Qk+l p Qk , 

a-k — Sk — Sk-i — Jk o,k Jk-i Ofc-i , 

Pk Pk-i 

or multiplying by Pk/a-k, and using (||), one gets the recurrence equation 

Pk = Qk+l fk — Qk — ^— fk-1 = Qk+l fk ~ rk fk-1 ■ (7) 

Pk-i a k 

Using @, Gosper gives an upper bound for the degree of / in terms of the degrees of p k , qk, 
and rk which yields a fast method to calculate fk, so that we finally find Sk, given by (||). This 
shows in particular that whenever ak possesses a closed form antidifference Sk then necessarily 
Sk is a rational multiple of a^: 

Qk+l fk 

s k = R k a k with R k = . 

Pk 

If any of the steps to find the polynomial fk fails, the algorithm proves that no hyper geometric 
term antidifference Sk of exists. 

Therefore the Gosper algorithm is a decision procedure which either returns "No closed form 
antidifference exists" or returns a closed form antidifference Sk of ak, provided one can decide 
the rationality of ak/a>k-l, i- e - one finds polynomials Uk,Vk such that a^/cLk-l = u k/ v k- in so 
far, the Gosper algorithm is an algorithm with input Uk and Vk rather than a^. 
Since without preprocessing, the user's input is ak rather than the polynomials Uk and Vk, the 
success of an implementation depends heavily on an algorithm quickly and safely calculating 
(uk,Vk) given a k - In Algorithm [j], we present such a method. It turns out that none of the 
existing implementations of Gosper's algorithm uses such a method, examples of which we 
will consider later. 
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In case, the Gosper algorithm provides us with an antidifference Sk of a^, any sum 

n 
k=m 

can be easily calculated by an evaluation of s/% at the boundary points like in the integration 
case. Note, however, that the sum 

n 

E(2) («) 

k=0 

e. g. is not of this type as the summand ( ^ J depends on the upper boundary point n explicitly. 
This is an example of a definite sum that we consider in § |3| 

Gosper implemented his algorithm in the MACSYMA nusum command, an implementation of 
the algorithm is distributed with the sum command of the Maple system (to check its use set 
inf olevel [sum] :=5), and one was delivered with Mathematica Version 1.2 (in the package 
Algebra/GosperSum.m). Another MATHEMATICA implementation was given by Paule and 
Schorn (1994). 

On the lines of Koornwinder (1993), together with Gregor Stolting I implemented the Gosper 
algorithm in Reduce (Koepf, 1994) and Maple, using the simple decision procedure for 
rationality of hyper geometric terms described in Algorithm [l] below rather than internal sim- 
plification procedures (like Maple's expand). In § ||, we will show that this makes our 
implementations that can be obtained from the author much stronger than the previous ones. 
It is almost trivial but decisive that the following is a decision procedure for the rationality of 
O'k/O'k-i for input (at least) of a special type: 

Algorithm 1 (simplif y_combinatorial) 

The following algorithm decides the rationality of a^jak-Y- 

1. Input: a\. as ratio of products of rational functions, exponentials, factorials, T function 
terms, binomial coefficients, and Pochhammer symbols that are rational-linear in their 
arguments. 

2. (togamma) 

Build afc/afc-i, and convert all occurrences of factorials, binomial coefficients, and Poch- 
hammer symbols to T function terms. 

3. (simplif y_gamma) 

Recursively rewrite this expression according to the rule 

r (a + k) = (a) k • T (a) 

((°)fc := a ( a + 1) • • • (a + A; — 1) denoting the Pochhammer symbol) whenever the argu- 
ments a and a + k of two representing T function terms have positive integer difference 
k. Reduce the final fraction cancelling common T terms. 

4. (simplif y_power) 

Recursively rewrite the last expression according to the rule 

b a+k = b k b a 

whenever the arguments a and a + k of two representing exponential terms have positive 
integer difference k. Reduce the final fraction cancelling common exponential terms. 
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5. The expression ak/ak-i is rational if and only if the resulting expression in step 4 is 
rational u^/vk, Uk,Vf. G Q[k]. 

6. Output: (uk,Vk)- 

Note that this result follows immediately from the given form of a& and therefore of the 

expression dk/ak-i considered. 

As an example, the rationality of a,k/a k -i of 

T(2k) 

ak - 



4 k T (Jfe)T (k + 1/2) 

is recognized by the given procedure, and from the resulting information {ak/dk-i = 1), by 
induction = 1/(2-^/7?) (Abramowitz and Stegun, 1964, (6.1.18)). 
Algorithm [j] does also apply to 

a fe = r(2£0-a4 fc r(£;)r(A; + l/2) , 

and leads to 

2 (2 k - 1) (k - 1) 



which is true whenever a ^ If ct = however, ak = 0, and therefore ak/ak-i is not 

properly defined. 

For this reason, Gosper's algorithm with input 

a k = T (2k) -J—A k T(k)T(k + 1/2) 

2^7T 

fails to find the true antidifference Sk = 1. 

We note, however, that in most cases also sums of ratios of the described form can be treated 
by the same method without using multiplication formulas of the T function explicitly. An 
important family of examples of this type will be considered in the next section. 
We note, finally, that an implementation of Algorithm [l] in a computer algebra system allows 
the user to enter his input in the form in which it is found in the literature, and no preprocessing 
is necessary. The user can be sure that the rationality is decided correctly. Unfortunately, 
this is not so with any of the current implementations. 

Maple's expand command which is used for that purpose in Zeilberger's (1990) and Koorn- 
winder's (1993) implementations of Zeilberger's algorithm (see § ^), e. g., does not the job 
required. The same is valid for Mathematical FactorialSimplif y procedure that comes 
in the package DiscreteMath/RSolve .m. Also, the Gosper implementation which comes with 
Maple's sum command, has the same failure. Paule and Schorn's Mathematica implemen- 
tation aborts in those cases with the message input not interpretable, whereas Gosper's 
nusum command gives the (wrong) error message errexpl NON-RATIONAL TERM RATIO TO 
NUSUM. 

Therefore in all these implementations, the fact that Gosper's algorithm is a decision proce- 
dure, unfortunately is completely lost. An example for that fact is given by the expression 



ak 



2 n V k J 



>n-l V k 



n — 1 
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with 

a k (n — k + 1) in — 2 k) 
a>k-i k (n — 2 k + 2) 

for which our Gosper implementations succeed very quickly, whereas Maple's sum command 
as well as Mathematical GosperSum, Paule/Schorn's Gosper, and Gosper's nusum fail if the 
input is not preprocessed by the user. 

This is a very simple example of an important type of examples considered next. 

3 The Wilf-Zeilberger Method 

Examples for an application of the Gosper algorithm in connection with Algorithm [l] are 
given by the Wilf-Zeilberger method on definite summation (Wilf and Zeilberger, 1990), see 
also Wilf (1993). 

The Wilf-Zeilberger method is a direct application of Gosper's algorithm to prove identities 
of the form 

s n := ]T F(n, k) = 1 (9) 

for which F(n, k) is a hypergeometric term with respect to both n and k, i. e. 

F(n,k) F(n,k) . , p 

and — — — r are rational functions with respect to both n and k, 



F(n-l,k) F(n,k-1) 

where n is assumed to be an integer, and the sum is to be taken over all integers k € 1L. 

To prove a statement of the form @ by the WZ method^], one applies Gosper's algorithm to 

the expression 

a k := F{n, k) — Fin — 1, k) 
with respect to the variable k. If successful, this generates G(n, k) with 

a k = F(n, k) - F(n - 1, k) = G(n, k) - G(n, k - 1) , (10) 

and summing over all k leads to 

s n - s n _! = £ (F(n, k) - F(n - 1, kj) = £ (G(n, k) - G(n, k - 1)) = 
kel. fceZ 

since the right hand side is telescoping. Therefore s n is constant, s n = s , and if we are able 

to prove so = 1, we are done. Note that so = 1 generally can be proved if the series considered 

is terminating in which case also no questions concering convergence arise. 

Since the WZ method only works if n is an integer, we can try to prove the statements of 

Bailey's list in Table |l] anyway only if one of the upper parameters of the hypergeometric 

series involved is a negative integer. The extension to the general case is over the capabilities 

of the methods of this article, and must be handled by other means. 

Note that the rationality of a^/cLk-i for the WZ method is decided by Algorithm [l] since 

a k F(n,k) - F(n - l,k) F(n,k) 1 F (n,k) 



a k ^ Fin, k — 1) — F(n —l,k—l) F(n,k-1) i F(n-i,fc-i) 

V ' I \ ' I \ ' F(n,fc-1) 



1 Note that Wilf and Zeilberger use forward differences rather that downward differences, whereas we decided 
to follow Gosper's original treatment. There is no theoretical difference between these two approaches, though. 
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Table 2: The WZ method 



Theorem 


n 


R(n,k) 


Vandermondc 


—a 


(h -4- k) (—n A- k) 
n (c + n — 1) 


Saalschutz 


n 


('/i _|_ M ( — n -L k) (n A- k) 

n (c + n — 1) (1 + a + 6 — c — n + fc) 


Kummer 


-b 


-1- A- ^1 I — n -A- k\ 

n (a + 2n) 


Dixon 


—c 


(n -A- k) ( n -A- k) (h -A- k) 

n (a — b + n) (a + 2 n) 


VVcilSOIl 

Whipple 


—c 


(n -4- h\ ( 7i -1- h\ (h -U 

(-1 + a + 6 + 2 n) (-2 n + 1 + k) (-2 n + k) 


Whipple 


—c 


(a + k) (a — 1 — k) (— n + £;) 
n(2-2n-e + fc)(l-2n-e + fc) 


Dougall 


n 


(2 a - 6 - c - <i + 2 n) (a + k) i-n + k) (b + k) (c + k) (d + k) 


n (a + 2 k) (a — b — c + n — d — k) (a — d + n) (a — c + n) (a — b + n) 


Dougall 


—e 


(a + k) {-n + k)(c + k) (d + k) 
n (a + 2k) (a — c + n) (a — d + n) 


Whipple 


—e 


(d + k) (-n + k) (a + k) 
n(a + 2k) (a — d + n) 


Bailey 


n 


(a 2 + 2a-wa+na + 2-2w-2kw + 2ka+2k + 2kn) (a + k) (-n + k) 
(—w + a + n) n (a + 2 k) (w + n — 1) 


Bailey 


n 


(-2 b - 2 b 2 + 2 nb + ab - 1 + n - k) (a + k) (-n + k) (b + k) 
nb (1 + 2 b - n + k) (a - 2 b + 2 n - 2) (a - b + n) 


Bailey 


n 


(2 b + ab + 1 - n + 2 kb + k) (b + k) (-n + k) (a + k) 
nb (a + 2 k) (1 + 2 b - n + k) (a - b + n) 



(a + k) (-n + k) (b + k) 
f n — — 

nb (a + 2 k) (2 + 2 b - n + k) (a - 2 b - 3 + 2 n) (a - b + n) 

■ (-8b-4b 2 + 6nb-ab-2n 2 +2nba-4 + 6n-2b 2 a + a 2 b-6k-8kb-4b 2 k+4kn+4kbn+2kba-2k 2 ) 
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Note moreover, that the application of Gosper's algorithm may be slow. But as soon as we 
have found the function G(n, k), we easily can calculate the rational function 

G(n,k) 

R{n,k) := 

F(n, k) 

by an application of Algorithm [|. R(n, k) is rational since the proof of Gosper's algorithm 
shows that G(n,k) is a rational multiple of a k = F(n,k) — F(n — l,k), G(n,k) = r(n,k) ■ 
(F(n, k) — F{n — 1, k)), say, so that 

itfn, «) = -=7 — — = r(n, fc) — - — — = r(n, fel l ^7 — — 

is rational. R(n,k) is called the rational certificate of F(n,k). Once the rational certificate 
of a hypergeometric expression F(n, k) is known, it is a matter of pure rational arithmetic 
(which is fast) to decide the validity of (||) since the only thing that one has to show is (10) 
which after division by F(n, k) is equivalent (modulo an application of Algorithm |l|) to the 
purely rational identity 

, , s F(n,k-l) F(n-l,k) , s 

l-R(n,k)+R(n,k-l)— ^ — — K — — = . (11) 

v ' ; v ' 1 F(n,k) F(n,k) v ' 

As an example, to prove the Binomial Theorem (compare (§)) in the form 

n n -I 

Sn :=E^^) = E F (I)= 1 (12) 

fc=0 fc=0 

by the WZ method, Algorithm [l] yields 

a k _ F(n,k) - F(n-l,k) _ (n - k + 1) (n - 2 k) 
a k -i ~ F(n, k - 1) - F(n - 1, Jfc - 1) ~ k(n-2k + 2) 

so that Gosper's algorithm can be applied, and results in 



G <"^^y( 2 (V) -(*))• 



This proves (|12| ) since sq = J2 1 = 1- 

k=0 

The rational certificate function is 



G(n,fc) fc-n 

/tin, k) = — — — = 

v ; F(n,k) n 



and the verification of identity (12) is therefore reduced to simplify the rational expression 

1-R(n,k)+R(n,k- i) ^"'*" 1 ) F (n~hk) = 1 fc-n | k-l-n k 2{n-k) 

F(n,k) F(n,k) n n n+l — k n 

to zero. 

Table § is a complete list of those identities of Bailey's list (Table |l]) that can be treated by 
the given method together with their rational certificates with which the reader may verify 
them easily. 
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Note that neither the statements of GauB and Bailey of argument x = 1/2 (p. 11) are accessible 
with respect to any of the parameters involved, nor can Watson's Theorem (p. 16) be proved 
by the WZ method with respect to Watson's original integer parameter a, nor can the method 
be applied to Whipple's Theorem (p. 16) concerning parameters a or 6 since in all these cases 
the term ratio ak/dk-i is not rational. 

Note further that both our Reduce and our Maple implementations generate the results of 
Table ||, and only the calculation of the rational certificate of Dougall's Theorem needs more 
than a few seconds. In the appendix, we will present some example results. On the other hand, 
Maple's sum command was not successful with a single example without preprocessing the 
input (entered in factorial or T notation), in most cases quickly responding with the incorrect 
statement Gosper's algorithm fails, and unsuccessfully trying other methods afterwards. 
In § ||, we consider a generalization of the WZ method. To be able to consider the most general 
case, we present an extended version of Gosper's algorithm next. 

4 An Extended Version of the Gosper Algorithm 

Here we deal with the question, given a nonnegative integer m, to find a sequence Sk for given 
afc satisfying 

Uk = Sk — Sk-m (13) 
in the particular case that Sk is an m-fold hypergeometric term, i. e. 

— — is a rational function with respect to k . (14) 

Sk~m 

Note that in the given case the input function itself is an m-fold hypergeometric term since 
by (H) and (0) 



1 



a k s k s k— m Sfc-m u k 



Oik—m Sk—m Sk—2m 1 s , ^k 

°K — m 

is rational, i. e., Uk and Vk can be chosen to be polynomials. 

Assume first, given a^, we have found Sk with Sk — Sk- m = ctk- Then we can easily construct 
an antidifference Sk of ak by 

Sk ■= Sk + s k -i-\ h s fc _ (m _i) (15) 

since then 



Sk — Sk-l — (Sfe H h S fc _( m _!)) — (Sfc_l H + Sk- m ) — Sk — Sfc-m — «fc • 

Note, however, that in general, this antidifference is not a hypergeometric term, but is a finite 
sum of hypergeometric terms. 

Assume next that given there exists a hypergeometric term with — Sfc- m = o-k- Then 
obviously also 

Sk Sk Sk-l Sfc-(m-l) 

Sk—m Sk—l Sk—2 Sk—m 

is rational, and therefore our algorithm below will find s&. 

An m-fold antidifference always can be constructed by an application of Gosper's original 
algorithm in the following way: 



10 



Algorithm 2 (extended_gosper) 

The following steps generate an m-fold antidifference: 

1. Input: afc, and m G IN. 

2. Define 6 fc := a km . 

3. Apply Gosper's algorithm to bk with respect to k. Get the antidifference i& of or 
the statement: "No hypergeometric term antidifference of and therefore no m-fold 
hypergeometric term antidifference of exists." 

4. The output s\~ := t k / m is a solution of fli~3| ) with the property (|l~4|). 
Proof: This is valid since a& = — Sk- m is equivalent to 

bk — O'km — Sfcm Skm—m — tk ^fc— 1 

and since 

tk Skm 



t-k—1 Skm—m 

describes the transformation between i& and Sk- □ 

As an example, we consider ak '■= k an d m = 2. Then bk = &2k = 2k k\, and Gosper's 

algorithm yields tk = 2{k + \)k\. Therefore Sk = tk/2 = (k + 2) ^|^! has the property that 

Sk — Sk-2 = ak ■ 
By (|i"5|), we moreover find the antidifference 

s k = s k + flfc-i = (k + 2) f ^ \ + (k + l) 



of a fe 

fc/3 



We consider two other examples: If = ( ' ^ then our algorithm generates the antidiffer- 



ence 



-*=37TXiT l(fc + 3)(>) + (fc + 2)( » ) + (* + r 



3(n + l 

and if«fc = ( fc / 2 )-( fc/2 n _ 1 ) then 



(2n + 3-k)(n + l-k) 
Sk 2(n + 2-k){n + l-k) 
(n + 2 - fc)(2n + 2 - k) 



f n \ ( n 

[ k=l I ~ [ fc-3 
\ 2 / V 



2 

n 



2(n + 2- k)(n + l-k) \\ k/2 J V fc/2 - 1 

Note, however, that we will use m-fold hypergeometric antidifferences rather than non-hyper- 
geometric antidifferences in the later chapters. 

Now, we give an algorithm that finds an appropriate nonnegative integer m for an arbitrary 
input function given as ratio of products of rational functions, exponentials, factorials, 
T function terms, binomial coefficients, and Pochhammer symbols that are rational-linear in 
their arguments: 
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Algorithm 3 (f ind_mf old) 

The following is an algorithm generating a successful choice for m for an application of Algo- 
rithm 

1. Input: a/; as ratio of products of rational functions, exponentials, factorials, T function 
terms, binomial coefficients, and Pochhammer symbols that are rational-linear in their 
arguments. 

2. Build the list of all arguments. They are of the form Pj/qj k + ay with integer pj and 
Qj, Pj/lj i n lowest terms, qj positive. 

3. Calculate m := lcmj^j}. 

Proof: It is clear that the procedure generates a representation for b k = a^ m with the given 
choice of m which is integer- linear in the arguments involved. Since in this case b^jb^-x is 
rational, Algorithm || is applicable. □ 

We mention that in our example cases above, the given procedure yields the desired values 

fc/3 



in 



2 for o fc := k (§)l, m = 3 for a k = ( n ) , and m = 2 for a fc = ( fc / 2 ) _ ( j~/2- 1 ) ' 



5 Extension of the WZ method 

In this section we will give an extended version of the WZ method which resolves the questions 
that remained open in § ^ so that finally Bailey's complete list (Table [l]) can be settled using 
a unifying approach. 

Assume that for a hyper geometric identity the WZ method fails. This may happen either 
because ak/a k -\ is not rational, or because there is no single formula for the result like in 
Andrews' statement 



3-^2 



— n , n + 3o , a 
3a/2 ,(3o + l)/2 



2 x , if n ^ (mod 3) 

l > = ' (nim3a + l) n ° therWiSe (16) 



which — together with many similar statements — can be found in a paper of Gessel and Stanton 
(1982), Equation (1.1). 

In such cases, we proceed as follows: To prove an identity of the form 

s n := F(n, k) = constant (n mod m constant) , (17) 
fceZ 

m denoting a certain positive integer, F(n, k) being an (m, l)-fold hypergeometric term with 
respect to (n, k), i. e. 

F(n,k) , F(n, k) . , p 

and — - — — are rational functions with respect to both n and k, 



F(n — m, k) F(n, k — I) 

and n assuming to be an integer, we apply our extended version of Gosper's algorithm to find 
an /-fold antidifference of the expression 

afc := F(n, k) — F(n — m, k) 
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Table 3: Gessel and Stanton's hypergeometric identities 

Eq. Identity 



(1.1) 
(1.2) 

(1.3) 
(1.4) 
(1.5) 
(1.6) 
(1.7) 



3^2 
5-^4 

3F2 
7-^6 



(1-8) 7 F 6 



(2d + 1) 2 „ (a + 6 + l/2)„ (a - b + 1)„ (a + b + l/2)„ (a - b + 1)„ (d + l/2)„ (d + 1)„ 



(3-7) 

(5.21) 
(5.22) 
(5.23) 

(5.24) 

(5.25) 
(5.27) 



2-fl 

3F2 
2F1 
2F1 

2F1 

2F1 
4F3 



—n , n + 3a , a 
3a/2 ,(3a+l)/2 



if n ^ (mod 3) 

n! (a + l) n/3 



(n/3)!(3a + l)„ 



otherwise 



2a ,26 ,1 - 2b ,1 + 2a/3 , -n 
a-6+l,a + 6+l/2, 2a/3 , 1 + 2a + 2n 

a, 6,a+l/2-6, l+2a/3,-n 
2a+l-26,26, 2a/3, l + a+n/2 



1 



(a + l/2)„ (a + l) n 
(a + 6+l/2)„(a-6+l)„ 



n! (o + l) n/2 2- n 



1/2 + 3a , 1/2- 3a ,-n 
1/2 ,-3n 



1 + 3a , 1 — 3a , — n 
3/2,-l-3n 

2a , 1 — a , — n 
2a + 2 ,-a- 1/2- 3n/2 



(|)!(a-6+l)„ /2 (6+±)„ /2 
3\ = (l/2-a) n (l/2 + a) n 
4 J (l/3)„ (2/3)„ 

_ (1 + a) n (1 - a) n 
(2/3)„ (4/3)„ 



if n odd 
otherwise 



1 



((n + 3)/2)»(n + l)(2a+l) 
(1 + (n + 2a + l)/2)„ (2a + n + 1) 



2a , 2b , 1 - 26 , 1 + 2a/3 , a + d + n + 1/2 , a - d , -n 
a - 6 + 1 ,a + 6+ 1/2 ,2a/3 , -2d - 2n , 2d + 1 , 1 + 2a + 2n 
a + 1) 2 „ (6 + d + l/2)„ (d - 6 + l) n _ (a + l/2)„ (a + 1)„ (6 + d + l/2)„ (d - 6 + l) r 



a, 6, a + 1/2 - 6,1 + 2a/3 , 1 - 2d , 2a + 2d + n , -n 
2a - 2b + 1 , 26 , 2a/3 , a + d + 1/2 , 1 - d - n/2 , 1 + a + n/2 



(6 + d) n/2 (d - 6 + a + l/2) n/2 n! (a + l) n/2 2"" 



(6 + l/2) n/2 (a + d + l/2) n/2 (d) n/2 (n/2)! (a - 6 + 1) 



if n odd 
otherwise 



-n ,-2n- 2/3 
4/3 

3a + 1/2 ,3a+l ,-n 
6a +1 ,-n/3 + 2a+l 



(5/6) n 
(3/2)„ 



/2 

(-27)' 



•1/2 



if n ^ (mod 3) 

(l/3) n/3 (2/3) n/3 



(1 + 2a)„ /3 (-2a)„/ 3 



otherwise 



(l/2)„ /27 V 



-n ,1/2 1 \ _ 

2n + 3/2 1) (2n + 3/2)„ \~I 

% - In \ 

8 = (-27)™ 



-n ,-1/3- 2n 
2/3 



-n ,n/2 + 1 
4/3 

-n ,1/2 
(n + 3)/2 



(l/2)„/2 



if n odd 

(-3)" ( ™ /2) otherwise 



(7/6) n/2 

if n odd 

(l/2) n/2 (3/2) n/2 



(5/6) n/2 (7/6) n/2 
1/3 - n , -n/2 , (1 - n)/2 , 22/21 - 3n/7 
5/6, 4/3, 1/21 -3n/7 
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otherwise 
27 



(-8)" 
1 -9n 



Bailey p. 



Table 4: The extended WZ method 

to R(n,k) 



11, Gaufi 



(b + k)(n- fc) 
{-b + n- l- 2k)n 



11, Bailey 



—a 



(2n-l)(n-fe) 
(c + n- l)(n + fc) 



16, Watson —a 



2 -2 



(c+fc) (6 + jfc) (n-fc) 



(-1 + n + 2 c) (-6 + n - 1 - 2 fc) n 



16, Whipple —a 



2 2 



(2n- l)(n-fc) (c+fc) 
(2 c - e + n) (-1 + n + e) (n + k) 



G.-S. Eq. 



TO 



R(n, k) 



1.1 



3 3 



(a + fc)(n-fc)(3a + 2n-3) 
(n + 3a + fc-2)(n + 3a + fc-l)n 



1.2 



(2 a + fc) (n - fc) (2 6 - 1 - k) (2 6 + fc) 
2n(3fc + 2a)(2a + 26-l + 2n)(a-6 + n) 



1.3 



2 4 



(6 + fc) (a + fc) (2 a - 2 & + 1 + 2 fc) (n - fc) 
n (3 fc + 2 a) (2 6 - 1 + n) (2 a - 2 6 + n) 



1.4 



1 3 



(n - fc) (6 a - 1 - 2 fc) (6 a + 1 + 2 fc) 
(12n-4fc)(3n-l-fc)(3n-2-fc) 



1.5 



1 3 



(n-fc) (3a-fc-l)(3a + l + fc) 
(3 n - 1 - fc) (3 n - fc) (3 n - k + 1) 



1.6 



(-4a+4an+18n 2 -20n+2-16nfc) (n-fc) (a-fc-1) (2a+fc) 
n(2a+l + 3n-2fc)(2a-l + 3n-2fc) (2 a-3+3 n-2 fc) (n- 1) 



1.7 



(2 a - 1 + 4 n + 2 d) (a - d + fc) (2 a + fc) (n - fc) (2 6 - 1 - fc) (2 6 + fc) 
n(2a+3fc) (2d+2n-fc) (2d+2n-l-fc) (2a+2 6-l + 2n) (a-b+n) 



1.8 



2 8 



(2rf-l-fc) (6+fc) (n-fc) (a+fc) (2 a-2 6+2 fc+1) (a + n+d- 1) 
n(2a+3fc) (-2 + 2d+n-2fc) (2 6-1+n) (2a-2 6+n) (2 a + 2 d+n+k- 1) 



3.7 



1 4 



(n - fc) (6 n + 2 - 3 fc) (7 n - 1 - 3 fc) 
(3n+l)(l + 2n)n 



5.21) 



3 2 



(3 a + 1 + fc) (6 a + 2 fc + 1) (n - fc) 
n (6 a + n) (-n + 6a + 3 + 3fc) 



5.22) 



(5 + 6fc)(l + 2fc) (n-fc) 
(24n + 4) (6n- l)n 



5.23) 



1 4 



(21 n - 7 - 9 fc) (6 n + 1 - 3 fc) (n - fc) 
(6n+ 1) (3n- l)n 



5.24) 



2 4- 



1 +3n 



5.25) 



(4n-4fc)(l + 2fc) (2 + 3fc) 
n(3n- 1) (1 + 3n) 



5.27) 



1 81 



(n - 1 - 2 fc) (n - 2 fc) (-1 + 3 n - 3 fc) 
n(3n- 1) (-l + 9n-21fc) 
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with respect to the variable k. (In most cases I = 1, so that Gosper's original algorithm is 
applied.) If successful, this generates G(n, k) with 



a k = F(n, k) - F(n - rn,k) = G{n, k) - G{n, k - I) , (18) 
and summing over all k leads to 

s n - s n . m = (Fin, k) - F(n - m, kj) = £ (G(n, k) - G(n, k - /)) = 
fceZ feeZ 

since the right hand side is telescoping. Therefore s n is constant for constant n mod m, and 
these constants can be calculated using suitable initial values. This can be accomplished if 
the series considered is terminating. Note, that again, the function 

acts as a rational certificate function. Once the rational certificate is known, it is a matter 
of pure rational arithmetic to decide the validity of (|l7]) since the only thing that one has to 
show is (18) which after division by F(n, k) is equivalent to the purely rational identity 

„. „, „.F(n,k — l) F(n — m,k) 

1 - R(n, k) + R(n, k - v ' - ; K — — -f^ = . 

v ' ; v 1 F{n,k) F(n,k) 



As an example, we prove (|lq) : In the given case, we set m := 3, I := 1, further 

_ (-n) k (n + 3a) k (a) k (n/3)! (3a + l) n ^3^ k 
[H,h fc! (3o/2) fc ((3o + l)/2) fc n!(o+l) n/3 

and notice that 

F(n,k) , F(n,fe) 

and 



F(n,fc-1) F(n-3,fc) 

are (complicated) rational functions (Algorithms || and ||). An application of Gosper's algo- 
rithm is successful, and leads to the rational certificate 

(a + k) (n - k) (3a + 2n - 3) 



R(n,k) 



(n + 3a + k-2)(n + 3a + k-l)n ' 
Therefore 

n 

22 F(n, k) = F{n, k) = constant (n mod 3 constant) , 
fceZ ^=0 

and statement (^) follows using three trivial initial values. 

Table || lists the hypergeometric identities of the Gessel-Stanton paper (note the misprint 
in Equation (1.4)), and Table |] contains their rational certificates (|l9[), calculated by our 
implementations, together with the certificates of Bailey's list (Table [j]) to which the WZ 
method did not apply. 

Note that in all cases considered, I = 1, so that the original Gosper algorithm is applied. 
Note, moreover, that Gessel and Stanton were not able to present proofs for their statements 
(6.2), (6.3), (6.5), and (6.6): Table [5] contains proofs. 
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Table 5: Gessel and Stanton's open problems 

Eq. Identity 



(6.2) 7 F 6 [ 



a+ 1/2 ,a ,b ,1 - b ,-n , (2a + l)/3 + n , a/2 + 1 
1/2 ,(2a - b + 3)/3 , (2a + b + 2)/3 , -3n , 2a + 1 + 3n ,a/2 
((2a + 2)/3)„ (2a/3 + 1)„ ((1 + b)/3)„ ((2 - b)/3)„ 



1 = 



(6.3) 
(6.5) 
(6.6) 



5*4 ( 
2^1 ( 
2*1 I 



((2a - b)/3 + 1)„ ((2a + 6 + 2)/3)„ (2/3) n (l/3)„ 

((2a + 2)/3)„(2a/3 + l)„ 



a + 1/2 , a , -n , (2a + l)/3 + n , a/2 + 1 
1/2 ,-3n ,2a+ 1 + 3n ,a/2 
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(2/3)„(l/3) r 



—n , — n + 1/4 
2n + 5/4 

— n , — n +1/4 
2n + 9/4 



(5/4). 



(2/3)„(13/12)„ V3 5 



2 6\n 



(9/4)s 



(4/3) n (17/12)„ V3 5 



2 6 \" 



Rational certificates 



Eq. m i?(n,fc) 



(a - 1 + 3 n) (a + k) (2 a + 2 fc + 1) (n - fc) (6 - 1 - k) (b + k) 



(6.2) 1 6 



(a + 2 fc) (3 n - k) (3 n - 1 - k) (3 n - 2 - fc) (2 a - b + 3 n) (2 a + 6 - 1 + 3 n ) 
(6 a - 6 + 18 n) (n - fc) (2 a + 2 fc + 1) (a + fe) 



(6.3) 1 



(a + 2 k) (3 n - fc) (3 n - 1 - k) (3 n - 2 - k) 



(52 n 2 - 13n-21 - 56 fc + 16 nfc - 32 fc 2 ) (n - k) (An- 1 -4k) 



(6.5) 1 



(108n- 27) (3n- 1) (1 + 12n)n 



(52 n 2 + 39 n - 55 - 84 fc + 16 nfc - 32 fc 2 ) (4 n - 1 - 4 fc) (n - fc) 



(6.6) 1 



(108 n- 27) (l + 3n) (5 + 12n)n 
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Similarly as the original WZ approach, our method is not capable, however, to prove Gessel- 
Stanton's (6.1), a non-terminating version of (6.2). Also, Gessel- Stanton's result (1.9) 



3-^2 



-sb + s + 1 , b — 1 , — n 
b + 1 , s(—n — b) — n 



(1 + s + sn) n b(n + 1) 
(l + s(6 + n)) n (b + n) 



is over the capabililites of our method since in this case the summand is an (m, /)-fold hyper- 
geometric term only for fixed (rational), but not for arbitrary s. 

Note, that our method not only unifies the proof of hyper geometric identities in a stronger 
fashion than the original WZ approach but moreover our Reduce and Maple implementa- 
tions do all the necessary computations completely automatically. We present some of the 
input and output in the appendix. 

Finally, we give examples of an application for which 1^1. To prove the identity (n G IN) 

-t(-*)* (;)•(?)-!. < 2 »> 

k=0 

we apply our extended WZ method with I = 2, and m = 1, and get the rational certificate 

_ (-fc + n-l)(-fe+n) 
K{n,k)- {n _ 1){ _ k + 2n _ 2) , 

which proves (^). Similarly one proves the statement (n G No) 

k=0 



6 The Zeilberger Algorithm 

In this section, we recall the celebrated Zeilberger algorithm (Zeilberger, 1990-1991), see also 
Graham, Knuth and Patashnik (1994) with which one can not only verify hyper geometric 
identities but moreover definite sums can be calculated if they represent hypergeometic terms. 
Zeilberger's algorithm deals with the question to determine a holonomic recurrence equation 

J 

£P»£(n-j) = (21) 
j=0 

with polynomials Pj in n, for sums 

S(n):= £F(n,fc) (22) 
fceZ 

for which F(n, k) is a hypergeometric term with respect to both n and k. 
Zeilberger's idea is to apply Gosper's algorithm in the following non-obvious way: Set 

J 

a k ■= F(n, k) &j («) F(n - j, k) 

3=1 
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with yet undetermined variables aj depending on n, but not depending on k. Then 

F(n,k)+ia,(n)F(n-j,k) 1 + £ * j (») 

a k j=i F(n,k) j=i ( ' 



Qk ~ l F(n,k-l)+j:^(n)F(n-j,k-l) l) l+S^ lff 

3=1 3=1 

turns out to be rational with respect to k, so the Gosper algorithm may be applied. 
If an application of Gosper's algorithm is successful it provides us with s k depending on n, 
and a set of rational functions <Jj{n) (the coefficients of f k are determined together with the 
unknowns cjj) such that 

J 

s k - = a k = F(n, k) + E Oj-(n) F(n - j, k) , 

3=1 

so that by summation 

E a k = E f F ( n > k ) + E °j W F ( n - ■?> fc ) 

fceZ fceZ \ 3=1 

J 

= S(n) + E s ( n _ j) = E ~ Sk ~i) = 
i=i fceZ 

since the right hand side is a telescoping sum. After multiplication with the common denom- 
inator this establishes the recurrence equation ( |2l| ) searched for. 

Koornwinder (1993) gives a rigorous description of Zeilberger's algorithm in the (most com- 
mon) case that the summation bounds are natural: a_i = a n+ \ = 0, i. e. the summation is 
for k = . . . n. 

Like for the Wilf-Zeilberger method, the Zeilberger algorithm is accompanied by a rational 
certification mechanism. 

Note that Zeilberger's algorithm can be applied to ratios of products of rational functions, 
exponentials, factorials, T function terms, binomial coefficients, and Pochhammer symbols 
that are integer-linear in their arguments with respect to both n and k. 

In the next section we will present a modified version of Zeilberger's algorithm that is appli- 
cable if the arguments of such expressions are rational-linear with respect to n and k. 
The application of Zeilberger's algorithm has the advantage over the WZ method that the 
right hand side of the hyper geometric identity does not have to be known in advance, but 
is generated by the algorithm (not to speak of the possibility to verify identities of other 
type). Therefore, Zeilberger's algorithm can be used to calculate definite sums rather than 
only verifying them. All identities mentioned in this article which could be verified with the 
WZ method, can be generated by Zeilberger's algorithm. 

Implementations of the Zeilberger algorithm were given by Zeilberger (1990) and Koorn- 
winder (1993) in Maple, and by Paule and Schorn (1994) in Mathematica. On the lines of 
Koornwinder (1993), we implemented the Zeilberger algorithm in Reduce (Koepf, 1994) and 
Maple, examples of which are given in the appendix. 
Note the following side conditions of the previous implementations: 

• Zeilberger: Here one must write the input into a file rather than on the command line. 
Supports only integer-linear input of a special form. 
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• Koornwinder: Supports only integer-linear input in hypergeometric notation. 

• Paule-Schorn: Supports only ratios of rational functions, products of exponentials, fac- 
torials, and binomial coefficients. 

Our implementations support the input in factorial-binomial-Gamma-Pochhammer as well as 
hypergeometric notation, and use Algorithm [l| for rationality decisions, and are therefore not 
bound to integer-linear input. 



7 An Extended Version of Zeilberger's Algorithm 

Our extended version of the Zeilberger algorithm deals with the question to determine a holo- 
nomic recurrence equation ( |2"l] ) for sums (^) for which F(n, k) is an (m, Z)-fold hypergeometric 
term with respect to (n, k), see § [|. 

In particular, this applies to all cases when the input function F(n, k) is given as a ratio of 
products of rational functions, exponentials, factorials, T function terms, binomial coefficients, 
and Pochhammer symbols that are rational-linear in their arguments with respect to both n, 
and k. 

First of all we mention that Zeilberger's algorithm may be applicable even though this is safely 
the case only if the arguments are integer-linear. An example of that type is the function 



E(n) := 



2-^1 



-n/2 , -n/2 + 1/2 
6 + 1/2 



A = g (-n/2) fc (-n/2 + l/2) fc 



fc!(6 + l/2) fc 

for which an application of Zeilberger's algorithm yields the recurrence equation 

(26 + n - 1) E(n) - 2(b + n - 1) E(n - 1) = , 
and therefore the explicit representation 

2" (b) n 



E(n) 



(26)„ 



Zeilberger's algorithm applies since F(n, k)/F(n — l,k) and -F(rz, k)/F(n, k — 1) are rational 
even though the representing expression for F(n, k) is not integer-linear in its arguments. 
On the other hand, not for every F(n, k) given with rational-linear T-arguments, the Zeilberger 
algorithm is applicable. An example for this situation is the left hand side of the Watson 
theorem with respect to variable a (see Table |]). 

We present now an algorithm which can be applied for arbitrary rational- linear input. 
Algorithm 4 (extended_sumrecursion) 

The following steps perform an algorithm to determine a holonomic recurrence equation ( |2l| ) 
for sums (f22|). 



1. Input: F(n, k), given as a ratio of products of rational functions, exponentials, factorials, 
r function terms, binomial coefficients, and Pochhammer symbols with rational-linear 
arguments in n and k. 

2. Build the list of all arguments. They are of the form Pj/qj n + Sj/tj k + ctj with integer 
Pj, qj, Sj,tj, Pj/qj and Sj/tj in lowest terms, qj and tj positive. 
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3. Calculate m := lcmj^j} and I := lcm{tj}. 

4. Define F(n,k) := F(mn,kl). Then F(n,k) is integer-linear in the arguments. 

5. Apply Zeilberger's algorithm to F(n,k). Get the recurrence equation 

J 

j=0 

with polynomials Pj in n, for the sum 

E(n) := ]T Hn,k) • 
feeZ 

6. The output is the recurrence equation 

J 

Pj(n/m) S(n — jm) = 

j=0 

for the sum 

S(n) := J] i^(n,fc) . 
fceZ 

Proof: Obviously our construction provides us with F(n,k) integer- linear in the arguments 
involved. Therefore Zeilberger's algorithm can be applied, and the result follows. □ 

As a first example, we apply our algorithm to the Watson function 

— n ,b ,c 



V ^ \ (-n + 6+l)/2,2c 



1 



with respect to the variable n to which Zeilberger's algorithm does not apply. In this case, 
the algorithm determines m = 2 and 1 = 1, and leads to the two- fold recurrence equation 

(b - 2c - n + 1) (n - 1) S(n - 2) - (6 - n + 1) (2c + n - 1) E(n) = 

from which the explicit right hand representation listed in Table || can be deduced for integer 
n since for positive values of n the Watson sum is finite, and therefore 

S(0) = 1 , 

and 

It turns out that our method is applicable to all identities considered in this paper to which 
Zeilberger's original approach does not apply. 

For example, we consider the three major identities of the paper of Gessel and Stanton (1982): 
The evaluation of (1.7) 



E(n) := 7 F ( 



2a ,2b ,1 - 2b ,1 - 2a/3 , a + d + n + 1/2 , a - d , -n 
a-b+1 ,a + b+l/2 , 2a/3 , -2d - 2n , 2d + 1 , 1 + 2a + 2n 
(a + l/2) w (a + l) n (b + d + l/2) w (d - b + l) n 
(a + 6 + l/2) n (a - 6 + l) n (d + l/2) n (d + l) n 



1 
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is found by a direct application of Zeilberger's algorithm with respect to n, leading to the 
equivalent recurrence equation 

= (2a + 26 + 2ra- l)(a - b + n)(2d + 2n - l)(d + n)£(n) 

+(2a + 2n - l)(a + n)(26 + 2d + 2n - 1)(6 - d - n)£(n - 1) . 

On the other hand, the evaluation of (1.8) 

/ a , b, a+ 1/2 -6,1 + 2a/3 , 1 - 2d , 2a + 2d + n , -n 
£(n) .- 7^6^_ 26 + lj2& ^ 2o ^ 3ja + d+1 ^ 2 ^ 1 _ ti _ ji ^ 2jl + a + n y 2 

if n odd 

(6 + d) n/2 (d - 6 + a + l/2) n/2 n! (a + l) n/2 



[ (6 + l/2) n/2 (a + d+ l/2) n/2 (d) n/2 (n/2)! (a - 6 + l) n/2 



otherwise 



cannot be handled with respect to n using Zeilberger's algorithm, but the extended version 
leads to the equivalent 2-fold recurrence equation 

= {n - 1 + 2d + 2a)(26 - n- 2a) (n - 1 + 26) (n - 2 + 2d)E(n) 

+(n - 1 + 2d - 26 + 2a) (n - 2 + 2d + 26)(2a + n)(n - l)S(n - 2) . 

A direct application of Zeilberger's algorithm is possible, however, with respect to the other 
variables involved (even with respect to a). 
Gessel-Stanton's open problem (6.2) 



EfrO 



7^6 



a + 1/2 , a , 6 , 1 - 6 , -n , (2a + l)/3 + n , a/2 + 1 
1/2 , (2a - 6 + 3)/3 , (2a + 6 + 2)/3 , -3n , 2a + 1 + 3n , a/2 



1 



((2a + 2)/3) n (2a/3 + 1)„ ((1 + 6)/3) n ((2 - 6)/3) n 



((2a - 6)/3 + l) n ((2a + 6 + 2)/3) n (2/3)„ (l/3) n ' 
again, can be solved directly with Zeilberger's algorithm leading to the recurrence equation 

= (2a + 6 + 3n-l)(2a-6 + 3n)(3n-l)(3n-2)S(n) 

+ (2a + 3n - 1) (2a + 3n) (6 + 3n - 2) (6 - 3n + 1) E(n - 1) . 



Finally, as an example with Z ^ 1, we consider (p0|), again. Our algorithm generates m = 1 
and 1 = 2, and the recurrence equations 



E(n) - Efn - 1) 



for 







and 2E(n) + E(n - 1) = 

=(»)=-(:)■(?)■ 



and 



respectively. 



8 Deduction of hypergeometric identities 

Finally, we mention that with a good implementation of Zeilberger's algorithm and our ex- 
tension at hand, it is easy to discover new identities. Just for fun, we realized the pattern in 
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Andrews' statement (16), and tried to generate similar ones: It turns out that 



3-^2 



and 



3^2 



— n , n + 2a , a 
2a/2 , (2a + l)/2 



— n , n + 4a , a 
4a/2 ,(4a + l)/2 



if n odd 

( -l)"/ 2 (l/2) n /2 

(1/2 + a) n/2 





(V2)„/2 



otherwise 

if n odd 
otherwise 



(1/2 + 2a) n/2 

Another example of a more deductive strategy is: Applying Zeilberger's algorithm to the 
general 2 iq polynomial 

S(n) := 2 Fn 



/ a , — n 








V & 





e. g., leads to the recurrence equation 

(b-l + n) S(rt) + (-2n + ra + xa - b + 2 - x) S(n - 1) - (x - 1) (n - 1) S(n - 2) = . 

It is therefore hypergeometric only if the coefficient of X(n — 2) is identical zero, i. e. if x = 1. 
This gives Vandermonde's identity. However, the coefficient of S(n — 1) can be made zero 
(equating coefficients), if we choose x = 2, and b = 2a, in which situation we get 

(n + 2a - 1) S(n) - (n - 1) S(n - 2) = . 

Therefore we have deduced the identity 

if n odd 
otherwise 



iF x \ 



a , — n 
2a 





(V2)n/2 



(l/2 + o) n/2 

We see that this method, to some extent, can be a substitute for the ingenuity of people 
like Dougall, Bailey, Andrews, Gessel or Stanton to find hypergeometric sums which can be 
represented by single hypergeometric terms. 

We finally give a strange example to demonstrate that our method can be of great help to 
find new identities. 

We try to find all hypergeometric functions of the form 



2F1 



a , — n 
n + b 



for which a, b and x are constants with respect to n, and for which a recurrence equation with 

only two terms S(n — j) is valid. 

The recurrence equation for S(n) turns out to be 

= — (x — l) 2 (n — 1) (n — 1 + b) (n — 2 + b) (xn + n — xa — x + bx) S(n — 2) 
+ (n - 1 + b) P(n, a, b, x) S(n - 1) 

+x (2 n + b — 1) (2 n + 6 — 2) (n — a — 1 + b) (xn + n — xa — 2x — 1 + 6x) S(n) , 

where P(n, a, 6, x) denotes a very complicated polynomial of degree 2 in n that does not have 
a rational factorization. To receive a recurrence equation for which only two terms E(n — j) 
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different from zero occur, we may set the coefficient lists with respect to n of any of the factors 
occurring zero, and try to solve for a, b and x. Note that since the resulting equations systems 
are polynomial systems, by Grobner bases methods these can be solved algorithmically. 
In our case, we receive either x = 1, i. e. the recurrence equation 

= -(n-l + 6)(2n-a + 6-l)(6-a + 2n-2)E(n-l) 
+ (2 n + b - 1) (2 n + b - 2) (n - a - 1 + b) S(n) , 

or we are led to the Kummer identity, i. e. to the values b = a + 1 and x = — 1 with the 
recurrence equation 

-2 (n + a) £(n - 1) + (2 n + a) £(n) = . 

The only exception occurs when we set the coefficient list with respect to n of the factor 
P(n, a, 6, x) zero, leading to the Kummer case again, and to the second solution set 

{a = 1/2, b = 3/2, x 2 - 6x + 1 = 0} . 

For x = 3 ± 2\/2, we have the recurrence equation 

-4 (2n- 1) (2n+ l)E(n-2) + (4n- 1) (4n+ l)S(n) = 

leading to the closed form representations 



2-^1 



1/2, -n 
n + 3/2 



3 + 2\/2 



otherwise 



and 



2^1 



1/2, -n 
n + 3/2 



3-2V2 



f 2(5/4) (w _ 1)/2 (7/4) (n _ 1)/2 /- 
— — — i v ^— (1-V2) if n odd 

5 ( n / y )(n-l)/2 (l<3/Oj(n-l)/2 

(3/4) ra/2 (5/4) ra/2 
(7/8) n/2 (9/8) n/2 

f 2(5/4) (n _ 1)/2 (7/4) (w _ 1)/2 

5(ll/8) (n _ 1)/2 (13/8) (n _ 1)/2 
(3/4) n/2 (5/4) n/2 



(1+V2) if n odd 



otherwise 



(7/8) n/2 (9/8) n/2 

in particular, for even n, the values at x = 3 + 2\/2 and x = 3 — 2a/2 are rational and equal: 



2-fl 



1/2 ,-2n 
2n + 3/2 



)2n 
fc=0 



( 



4n + 2k 
2k 



) 



(3 ± 2^2) 



fc _ (3/4) ra (5/4) ra 
(7/8) n (9/8) n 



We will discuss the given method in greater detail in a forthcoming paper. 
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Appendix 

In this appendix, we give a short description of a Maple implementation, which I implemented 
together with Gregor Stolting on the lines of Koornwinder (1993), incorporating Gosper's and 
Zeilberger's algorithms and the extensions of this article, and present some of its results. Our 
Reduce implementation is described elsewhere (Koepf, 1994). 
After loading our package, one can use the following Maple functions: 
gosper(f ,k) determines a closed form antidifference. If it does not return a closed form 
solution, then a closed form solution does not exist. 



gosper (f ,k,m,n) determines 



J2 a * 

k=m 



using Gosper's algorithm. This is only successful if Gosper's algorithm applies. 

extended_gosper (f ,k,m) determines an m-fold hyper geometric antidifference. If it does 
not return a solution, then such a solution does not exist. 

sumrecursion(f ,k,n) determines a holonomic recurrence equation for 

oo 

summ(n) = f( n ^k) with respect to n if f(n,k) is hypergeometric with respect to 

k=— oo 

both n and k. The resulting expression equals zero. 
sumrecursion(f ,k,n, j) searches only for a holonomic recurrence equation of order j. 
extended_sumrecursion(f ,k,n,m, 1) determines a holonomic recurrence equation for 

oo 

summ(n) = J2 f{ n -> with respect to n if f(n, k) is an (to, Z)-fold hypergeometric term 

k=— oo 

with respect to (n,k). 



hyperrecursion (upper, lower, x,n) determines a holonomic recurrence equation with re- 
al, 0,2, a p 
v h, h, ■■■ , b q 

upper parameters, and lower = {6 1; b 2 , ... , b q } is the list of lower parameters depending 
on n. 



spect to n for p F q 



where upper = {ai, a-2, ■ ■ ■ , a p } is the list of 



hyperrecursion (upper , lower, x,n, j) searches only for a holonomic recurrence equation of 
order j. 

hyperterm (upper, lower, x,k) yields the hypergeometric term 

{oi)k ■ {02)k ■ ■ ■ {Op)k x k 



(h)k- {h)k---(b q ) k k\ 



with upper parameters upper = {ai, ■ ■ ■ , a p }, and lower parameters lower = 
{h,b 2 , ...,b q } 

simplif y_gamma(f ) simplifies an expression f involving only rational functions, exponentials 
and r function terms according to a recursive application of the simplification rule 
r (a + 1) = aT (a) to the expression tree, see Algorithm |l[ 
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simplif y_combinatorial (f ) simplifies an expression f involving exponentials, factorials, T 
function terms, binomial coefficients, and Pochhammer symbols by converting factorials, 
binomial coefficients, and Pochhammer symbols into T function terms, and applying 
simplif y_gamma and simplif y_power to its result. If the output is not rational, it is 
given in terms of T functions, see Algorithm [l|. 
The Maple function 

WZ : =proc ( summand , k , n , m) 
local tmp,gos; 

tmp : =summand-subs (n=n-m , summand) ; 

gos : =extended_gosper (tmp , k , m) ; 

RETURN (simplif y_combinatorial (gos/ summand) ) 
end: 

therefore, calculates the (m,Z)-fold rational certificate (|lG|) of F(n,k). 
Here are some results of the implementation: 

|\~/| Maple V Release 3 (FU-Berlin) 

._l\l l/l_. Copyright (c) 1981-1994 by Waterloo Maple Software and the 
\ MAPLE / University of Waterloo. All rights reserved. Maple and Maple V 

< > are registered trademarks of Waterloo Maple Software. 

I TyP e ? f° r help. 

> read summation; 

> # see (SIAM Review, 1994, Problem 94-2) 

> gosper((-l)~(k+l)*(4*k+l)*(2*k) ! / (k ! *4~k* (2*k-l) *(k+l) ! ) ,k) ; 

(k + 1) 

(2 k) ! (-1) 



k 

(k + 1) ! 4 k! 

> # Dougall 

> WZ (hyperterm({a , l+a/2 ,b , c , d , l+2*a-b-c-d+n , -n} , 

{a/2 , 1+a-b , 1+a-c , 1+a-d , 1+a- ( l+2*a-b-c-d+n) , 1+a+n} , 1 , k) / 
hyperterm({l+a, 1+a-b-c, 1+a-b-d, 1+a-c-d, 1}, 
{1+a-b, 1+a-c , 1+a-d, 1+a-b-c-d} ,l,n) ,k,n,l) ; 

(2 a - d + 2 n - c - b) (a + k) (- k + n) (b + k) (c + k) (d + k) 

n (a + 2 k) (a - b - c - d + n - k) (a - b + n) (a - d + n) (a - c + n) 

> sumrecursion(hyperterm({a, l+a/2 ,b, c ,d, l+2*a-b-c-d+n, -n} , 
{a/2 , 1+a-b , 1+a-c , 1+a-d , 1+a- ( l+2*a-b-c-d+n) , 1+a+n} , 1 , k) / 
hyperterm({l+a, 1+a-b-c, 1+a-b-d, 1+a-c-d, 1}, 

{1+a-b, 1+a-c, 1+a-d, 1+a-b-c-d}, l,n) ,k,n) ; 

summ(n) - summ(n - 1) 

> hyperrecursion({a, l+a/2, b,c,d, l+2*a-b-c-d+n, -n} , 

{a/2, 1+a-b, 1+a-c, 1+a-d, 1+a- (l+2*a-b-c-d+n) , 1+a+n}, l,n) ; 
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- (a + n) (a - c - d + n) (a - b - d + n) (a - b - c + n) summ(n - 1) 



+ summ(n) (a - d + n) (a - c + n) (a - b + n) (a - b - c - d + n) 

> # Gessel-Stanton (6.2) 

> WZ(hyperterm({a+l/2,a,b, l-b,-n, (2*a+l) /3+n, a/2+1} , 
{1/2 , (2*a-b+3) /3 , (2*a+b+2) /3 , -3*n, 2*a+l+3*n, a/2} , 1 ,k) / 
hyperterm ({ (2*a+2) /3 , 2*a/3+l , ( 1+b) /3 , (2-b) /3 , 1} , 
{(2*a-b)/3+l, (2*a+b+2)/3,2/3,l/3},l,n) ,k,n) ; 

6 (a - 1 + 3 n) (a + k) (2 a + 2 k + 1) (- k + n) (b - 1 - k) (b + k)/( 

(a + 2 k) (3 n - k) (3 n - 1 - k) (3 n - 2 - k) (2 a - b + 3 n) 

(2a+b-l+3n)) 

> sumrecursion(hyperterm({a+l/2,a,b, l-b,-n, (2*a+l) /3+n, a/2+1} , 
{1/2 , (2*a-b+3) /3 , (2*a+b+2) /3 , -3*n, 2*a+l+3*n, a/2} , 1 ,k) / 
hyperterm ({ (2*a+2) /3 , 2*a/3+l , ( 1+b) /3 , (2-b) /3 , 1} , 
{(2*a-b)/3+l, (2*a+b+2)/3, 2/3,1/3}, l,n) ,k,n) ; 

summ(n) - summ(n - 1) 

> hyperrecursion({a+l/2,a,b, l-b,-n, (2*a+l)/3+n, a/2+1}, 
{1/2, (2*a-b+3)/3, (2*a+b+2) /3, -3*n,2*a+l+3*n, a/2} , 1 ,n) ; 

- (3 n - 2 + b) (3 n - 1 - b) (2 a + 3 n) (2 a - 1 + 3 n) summ(n - 1) 

+ summ(n) (3 n - 1) (3 n - 2) (2 a - b + 3 n) (2 a + b - 1 + 3 n) 

> # The following two sums are identified to be equal, see Strehl (1993) 

> sumrecursion(binomial(n,k)~3,k,n) ; 

2 2 2 

- 8 (n - 1) summ(n - 2) - (7 n - 7 n + 2) summ(n - 1) + summ(n) n 

> sumrecursion (binomial (n,k) ~2*binomial(2*k,n) ,k,n) ; 

2 2 2 

- 8 (n - 1) summ(n - 2) - (7 n - 7 n + 2) summ(n - 1) + summ(n) n 

> simplif y_combinatorial ( (binomial (n , k) -binomial (n-2 , k) ) / 
(binomial (n-3 , k) -binomial (n-6 , k) ) ) ; 

(n - 5) (n - 4) (n-3) (n-2) (- k + 2 n - 1) 



2 2 
(3 n - 24 n - 3 k n + 12 k + k + 47) (n - 2 - k) (- k + n - 1) (- k + n) 

> extended_gosper (binomial (k/2,n) ,k,2) ; 
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(1/2 k + 1) binomial(l/2 k, n) 
n + 1 

> WZ(binomial(n,k)*binomial(k/2,n)*(-l)~k*(-2)~n,k,n,l,2) ; 

(- k + n - 1) (- k + n) 
(n - 1) (- k + 2 n - 2) 

> extended_sumrecursion (binomial (n,k) *binomial(k/2 ,n) *(-l)~k*(-2)~n,k,n,l,2) ; 

summ(n) - summ(n - 1) 

> extended_sumrecursion (binomial (n,k)*binomial(k/2,n)* (-2) ~n,k,n, 1 ,2) ; 

summ(n) - summ(n - 1) 

> hyperrecursion({-n,n+2*a,a},{2/2*a, (2*a+l)/2},2/4,n) ; 

(n - 1) summ(n - 2) + (n + 2 a - 1) summ(n) 

> hyperrecursion({-n,n+4*a,a},{4/2*a, (4*a+l)/2},4/4,n) ; 

- (n - 1) summ(n - 2) + (n + 4 a - 1) summ(n) 
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